The nonlinear electromigration of analytes into confined spaces 
[electromigration into confined spaces] 
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We consider the problem of electromigration of a sample ion (analyte) within a 
uniform background electrolyte when the confining channel undergoes a sudden 
contraction. One example of such a situation arises in microfluidics in the electroki- 
netic injection of the analyte into a micro-capillary from a reservoir of much larger 
size. Here the sample concentration propagates as a wave driven by the electric 
field. The dynamics is governed by the Nerst-Planck-Poisson system of equations 
for ionic transport. A reduced one dimensional nonlinear equation describing the 
evolution of the sample concentration is derived. We integrate this equation numer- 
ically to obtain the evolution of the wave shape and determine how the the injected 
mass depends on the sample concentration in the reservoir. It is shown that due to 
the nonlinear coupling of the ionic concentrations and the electric field, the concen- 
tration of the injected sample could be substantially less than the concentration of 
the sample in the reservoir. 

Keywords: electromigration dispersion, electrokinetic injection, nonlinear 

waves 



1. Introduction 

In microfluidic systems, electrokinetic injection is often e mployed i n order to insert 
the s ample from a reservoir into a microfluidic channel ([Landers! . Il99rl ICamillerl 



1998). A simplified problem depicting the electrokinetic injection process is shown 
in Figure [T] The reservoir initially contains the sample dissolved in a background 
electrolyte of positive and negative ions. The channel contains only the background 
electrolyte. At time t = 0, a voltage, V is applied to an electrode in the reservoir 
located far from the inlet, while the channel outlet (assumed infinitely far away) is 
electrically grounded. We would like to describe how the concentration of the sample 
evolves with time as the sample moves into the micro-channel. For simplicity, we 
will assume that the capillary walls are uncharged so that there is no electroosmotic 
flow. 

When the concentration of sample in the reservoir is low, the dynamics is linear 
and is described by an advection diffusion equation with constant advection speed. 
At high concentrations, the evolution is nonlinear due to the phenomenon of electro- 
migration dis persion. The physical m echanism of electromigration dispersion may 



be explained (jGhosal fc Chenl . l2010al) in the following way: when the concentration 
of sample ions is significant in comparison to that of the background electrolyte, 
the electrical conductivity of the solution is locally altered. However, charge con- 
servation and local electro neutrality requires the electric current to be the same 
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at all points along the axis of the capillary. Therefore, by Ohm's law, (ignoring, 
for the moment, the diffusive contributions to the current), the electric field must 
change axially, since, the product of the conductivity and electric field must remain 
constant. This varying electric field alters the effective migration speed of the sam- 
ple ions, which in turn, alters its concentration distribution. The transport problem 
for sample ions then becomes nonli near and shock li ke structures similar to those 
familiar from breaking water waves (|Whithaml . 1 19T4T ) can arise. 

The problem has a certain similarity with the classical "dam break problem" 
in hydrodynamics where a large stationary volume of water held back by a wall is 
sudd e nly released by partia l or co mplete removal of the restraining wall ( Whitham, 
19551 IStokerL Il948t iRitter . 1892 ). However, in the present problem, the driving 



mechanism is the variation of the electric field along the channel instead of hydro- 
static pressure variations due to changes in wave height. 

In the next section we present the mathematical formulation for a minimal 
problem consisting of three ionic species - the analyte, a co-ion and a counter-ion, 
where the diffusivities of the three species are equal, the charge of the analyte 
is constant, and, the background electrolyte is fully dissociated. We will see that 
even in this highly idealized limit, the nonlinearities inherent in the system lead 
to surprising behavior. The idealizations permit reduction of the problem to the 
solution of a single partial differential equation in one dimension. In section [3j this 
equation is integrated numerically, and observations are presented for the behavior 
of the sample ion concentration as it moves into the micro-channel. We find that at 
high concentrations, the concentration of sample ions that enter the micro-channel is 
significantly lower than that in the reservoir. In fact, as the concentration of sample 
ions in the reservoir is increased, the concentration within the injected plug in the 
channel does not increase indefinitely; but, approaches a limiting value determined 
by the valence of the analyte ion and background electrolytes. In section dj we 
identify the reason for this behavior and use simple arguments based on conservation 
laws to deduce the limiting value of the sample concentration in the channel. This 
is compared to results from the numerical simulations and good agreement is found. 
Our findings are summarized in section [5] and its implications in a broader context 
as well as its limitations are discussed. 



2. Problem Formulation 

We consider a three ion system consisting of sample ions, co-ions and counter- 
ions all of equal diffusivity (D). It follows, from the Einstein relation, that the 
mobility (u) of the three ionic species are also the same. However, the three species 
have different valences: z p ,z n ,z and therefore different electrophoretic mobilities: 
fip = z p eu, fi n = z n eu,fi — zeu, e being the proton charge. In this paper, the suffix 
p will generally indicate the positive ion (cation), n will indicate the negative ion 
(anion) and the absence of a suffix will indicate the sample ion (analyte). Thus, z p 
is positive, z n is negative and z could be of eith er sign. The discussion in the rest 
of this section follows closely the earlier work of Ghosal &: Chenl ( 2010af) . 



We consider the simplest situation where the ions migrate from a semi-infinite 
reservoir to a uniform channel with no wall charge. Then the problem is entirely 
one dimensional, and, the coupled equations describing the concentrations of the 
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Zone I 
(Reservoir) 



Cp— Cp* 



Zone II 
(Diffusion) 



Zone III 
(Advection) 



Zone IV 
(Background) 



Cp— Cp 

Cn = Cn° < 

c=0 



Kim = Kiv 



Figure 1. Schematic diagram illustrating the electrokinetic injection of a sample 
(concentration c* in the reservoir) into a micro-channel 
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(2.1) 
(2.2) 
(2.3) 



where E is the local electric field, x is the distance along the capillary and t is 
time. We will assume that the x-axis points in the direction of front motion from 
the reservoir into the capillary. Since characteristic spatial scales are always muc h 
larger than the Debye length, local electro-neutrality holds ( Ghosal fc Chenl . 2010h . 
Thus, 

z p Cp + z„Cn + zc = 0. (2.4) 

If we multiply equations (|2.1[) - (l2.3p by the respective ionic charges, sum them, 
and use equation (|2.4p , we get an equation that describes the constancy of electric 
current 

d 

— [e 2 u(z 2 c p + z 2 lCn + z 2 c)E] = 0. (2.5) 

Note, that, the net contribution from the diffusive fluxes vanish exactly on account 
of the assumption of equal diffusivity of ions and local electro-neutrality. Equation 
(|2.5p may then be integrated to yield 



(,ZpCp -\- z n c n 4- z c)E — (^Cp 4- z n c n )-£/oq, 



(2.6) 



where Eoo is the electric field in the capillary very far away from the inlet, and, 
, are respectively the corresponding negative and positive ion concentrations 



Article submitted to Royal Society 



The nonlinear electromigration of analytes into confined spaces 



■5 



in the background electrolyte. Clearly, and c£° are not independent but are 
related by the electro-neutrality condition 

z pCp ° + z n c™ = 0. (2.7) 

If we now introduce the Kohlrausch function (iKohlrauschL Il897h 

K(x,t) = (c p + c n + c)/u, (2.8) 

then it follows from equations (I2.1j) - (|2.3[) and (|2.4I) that K satisfies the diffusion 
equation 

The three algebraic relations: (12.41) . (I2.6[) and (I2.8[) . may be used to express all 
four dependent variables in the problem: c p ,c n ,c and E in terms of any one of 
them and the function K(x,t). We choose to express all the variables in terms of 
c and K. These relations take a particularly simple form if we use the normalized 
concentration 

Hi 

and 



(2.10) 



e ( x ^> = " : i^ - ( 2 - n ) 



zip 

In terms of these variables, we then have 



E = -j-^L-, (2.12) 



Z Z " ^ ^0, (2.13) 



(2.14) 



so that, equation (|2.3p gives the following evolution equation for 



3 ' "* ^-i>S, ("5) 



9f dx \9 — a<j> J dx 2 ' 

where v = zeuEoa is the velocity of an isolate d sample ion in the cons tant electric 
field E = Eoo, and, following the notation of iGhosal fc Chenl (|2010al) . a = \{z — 
z n )(z — Zp)/z n (z p — z n )], is the "velocity-slope parameter". 

Since 9 is proportional to K, it is also a passive scalar. The initial and boundary 
conditions on 9 in the domain < x < oo are 

9(x,0) = 1 ifx>0, (2.16) 
0(0, t) = 9*, (2.17) 
0(oo,t) = 1, (2.18) 

where 0* is the value of 9 in the reservoir. The reservoir is considered to be "well 
mixed" so that all ionic concentrations remain constant within it. An exact solution 
for 9 may then be found using Fourier's method, 



= 6* + (l-e*)eri[-=), (2.19) 
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where 

erf(x) = -= / e~ y dy (2.20) 
V 71 " Jo 

is the error function. Since the concentration of sample in the reservoir is being 
held constant at some value <fi = <fi* , we have 

g* = <_ + ZE^*_P (221) 

The superscript * indicates the value of the corresponding variable in the reservoir. 
Thus, c* is the concentration of negative ions in the reservoir. 

Equation (|2.15l) . on the other hand, does not readily admit an analytical solu- 
tion, since, unlike the equation for 9, it cannot be reduced to an ordinary differential 
equation by means of the similarity variable rj = x/V ADt. We will therefore obtain 
the time evolution of the normalized sample concentration (f>(x, t) by substituting 
equation (12.191) in equation (I2.15|) and integrating the resulting one dimensional 
partial differential equation numerically. Before we present this solution in the next 
section, it is useful to note that at large distances from the reservoir, or more pre- 
cisely, if x ^> %/ ADt, equation (|2.19l) implies that 9 ~ 1, so that equation (|2.15p 
reduces to 

~6t + (r=t^) = D d^ (2 - 22) 

which was studied earlier bv lGhosal fc Chenl (|2010al) . 



3. Numerical Simulations 

We will find it convenient to express our results in terms of a characteristic length, 
w (the channel width), which gives a characteristic time w/v. Then clearly, the only 
parameters in the problem are Pe = vw/D, which may be regarded as a "Peclet 
number" based on the electromigration speed, v, or a "field strength parameter"; 
the two valence ratios z n /z, Z p /z; and, the parameter 9* characterizing the ionic 
composition of the mixture in the reservoir. For definiteness, we will assume that 
the sample is a cation (z > 0), and that 

4 = (3.1) 

This will be true, for example, if the electrolyte in the reservoir is prepared by adding 
the sample in solid form (so there is no dilution) to a stock solution of the back- 
ground electrolyte that fills the micro-channel (since the sample contributes only 
counter-ions, the co-ion concentration is not altered). Using the electro-neutrality 
conditions, 

fl» = l- Zp ( Z ~_ Zn \ <i>*. (3.2) 

Zji \ Zp Z n J 



The initial and boundary conditions for <j) are 



4>{x,{)) = 0, if a; > (3.3) 
0(O,i) = 0*, (3.4) 
<p(oo,t) = 0. (3.5) 
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Time (vt/w) 



Figure 2. The increase of sample ions in the micro-channnel as a function of dimensionless 
time at different dimensionless reservoir concentrations, <f>* . The solid line corresponds to 
Pe = 10 and the dashed line corresponds to Pe = 200. 



We keep the valence ratios fixed at z p /z = 2, z n /z = — 1. 

Since we have reduced our problem to a one dimensional one, the numerical 
integration is quite straightforward. We used a finite volume method to discretize 
equation (I2.15[) in space using an adaptive grid refi nement algorithm tha t is en - 
abled by applying the Matlab library "MatMOL" (jVande Wouwer et all . 12004) . 



The spatially discretize d system of equations is the n integrated in time using the 



Matlab solver "ode45" (jShampine fe Reichehi 119971 ) which is based on an explicit 



Runge-Kutta (4,5) formula. 

Figure [5] shows the amount of sample that has moved into the reservoir at time 
t, defined as 



i r 

m(t) = - / 
w Jo 



(x, t) dx. (3-6) 



Simulations are shown for Pe = 10 and 200. At long times, or more precisely, if 
vt/uu 3> Pe^ 1 , the expected asymptotic form for m(t) is 

m(t) ~ {v4>*/w)t. (3.7) 

Figure [2] shows that indeed, at large times, m(t) oc t. We determined the slope 
of the m(t) curve at large times and plotted the dimensionless quantity wrh/v 
as a function of (j>* in Figure [3] It is seen that wrn /v ~ <fi* as long as <f>* -C 1 
but as the concentration of sample in the reservoir is increased, the injection rate 
does not increase indefinitely but tends to saturate to a limiting value, that is, 
wrh/v ~ (f>* m < 4>* . 

The reason for this saturation of the injection rate may be understood by ex- 
amining the concentration profiles <p(x, t) shown in Figure 01 It is seen, that as the 
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Figure 3. The normalized concentration of analyte in the micro-channel (<j>m) as a function 
of the corresponding quantity in the reservoir (<)>*). The solid line indicates the result 
expected from linear theory, the dashed line is the result predicted by the nonlinear theory, 
equation (|4.7|l . The symbols (circle: Pe=10, square: Pe=200) are obtained by numerical 
integration of equation (|2.15p . 



sample moves into the capillary, the concentration in the capillary is approximately 
equal to the the concentration in the reservoir as long as this concentration remains 
small (upper panel, <p* — 0.01). However, if the concentration in the reservoir is 
large (lower panel, (j>* = 1.0) the concentration in the channel approaches a limiting 
value «b* m independent of the reservoir concentration. 

In Figure [51 we have reproduced the concentration distribution <p(x,t) that ap- 
pears in the lower panel of Figure SI together with the concentrations of the other 
ionic species 4> p (x, t) and 4> n {x, t), the local electric field E{x, t), and, the local elec- 
trical conductivity due to ions, <j{x 1 i). These quantities are obtained from equations 
(|2.12p . (I2.13P and (|2.14l) . The mechanism for the reduction of the concentration from 
<j> = cf>* in the reservoir to <f> = <j)* n < (jf in the capillary may be understood from 
these graphs. At large values of the sample concentrations in the reservoir, there 
is a a sharp rise in the solution conductivity in the immediate vicinity of the inlet 
resulting in a drop in the electric field. The flux of sample into the micro-channel is 
then determined by a combination of this greatly reduced electromigrational flux 
and the diffusive flux and this is less than the expected flux indicated by equation 
(|3.7p . which should hold in the linear regime. 



4. Analysis 

An approximate theoretical determination of the critical concentration (\>* m may be 
provided using the conservation equations. The method of doing this i s in fact en- 
tirely analogous to that of the "Moving Boundary Equations" (MBE) ( Dole! . 1945f) 
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Figure 4. Concentration profiles in the micro-capillary at successive times: 
vt/w — 2,4,6,8, 10. If the reservoir concentration is large (lower panel), a concentration 
jump develops at the entrance of the capillary, with a narrow diffusion zone connecting 
the values in the reservoir with the much smaller value of the concentration within the 
capillary. The concentration jump is insignificant if the reservoir concentration is small 
(upper panel). The dashed line is <\>* m evaluated using equation l|4.8[) . Here Pe=200. 



for describing advancing fronts (e.g. in isotachophoresis). The conceptual framework 
is illustrated in Figure [T] The domain is decomposed into four parts: (I) the "Reser- 
voir Zone" within the well mixed body of the reservoir where all concentrations are 
constant (II) the "Diffusion Zone" which is the initial part of the micro-channel 
of length of order \j2Dt - here the ionic diffusive flux and the electromigrational 
flux driven by the electric field are comparable, (III) the "Advection Zone" where 
diffusive fluxes are negligible (because concentration gradients are small) and ionic 
fluxes are primarily due to electromigration, and, (IV) the "Background Zone" far 
downstream where the sample ions have not yet penetrated. 

The boundary between the Zones (I) and (II) is stationary (and located at the 
channel inlet) whereas the boundaries between the other zones propagate to the 
right. The arrows in Figure [T] indicate fluxes of ions across the stationary zone 
boundary. Equality of the fluxes for each species across the zone boundary require: 



7rr0* E {1) (r) 



(4.1) 
(4.2) 



where E^{r) represents the electric field in the reservoir at a radial distance r 
from the inlet. The superscript (I), (II), (III) or (IV) indicate the zone in which the 
corresponding variable is evaluated. Clearly qt>^ = 4>*, <j>n = <j)* n and cj^ 111 ^ = <p* m . 
Therefore, taking the ratio, 



(4.3) 



Article submitted to Royal Society 



10 



Z. Chen and S.Ghosal 




Figure 5. Same as in Figure [4] except that profiles of all dependent variables are shown but 
at a fixed time instant, vt/w = 2. The dashed lines are the concentrations in the capillary 
calculated using equation (|4.7[) . The bottom panel shows the electric conductivity (a) 
normalized by its downstream value (o~o). 

The Kohlrausch function, or, equivalently, 6(x,t), propagates only by diffusion. 
Therefore, 8(x,t) differs from unity only in Zones (I) and (II). Thus, 

^)+^/) +c = 1 (44) 

The electro-neutrality condition (valid in all zones) for Zone (III) is: 

^f'+zi'"'!*^. (4.5) 

By combining equations (|4.3p and (|4.4[) and using the electro-neutrality condition, 
equation (|4.5[) . we get an equation for determining cf>* n 

(1 - z/zp)^ + r(l - z n /z p )4>* m = l-z n /z p , (4.6) 

where the ratio <Pnl 't'* — r is a constant determined by the ionic composition of the 
electrolyte in the reservoir. Solving the above linear equation for <\>* m we have 

<t>*m = [r+ {z p - z)/(z p - z„)] _1 . (4.7) 
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In our numerical experiment we assume that the sample is a cation, and that 
the co-ion concentration in the reservoir is the same as in the channel, so that 0* = 
0£° = —z n /z p . The electro-neutrality condition then implies that 0* = l — (z/z n )</>*, 
so that equation (|4.7I) becomes 



where 



l + £0* 



Z<p{z Z-n) 



(4.8) 



Z-n { Zp Z n j 

is a positive number. An expression for 9* in terms of the sample ion concentration, 
0* , may also be derived by simple algebra: 

e* = i- Zp j z ~_ Zn \ <i>*. (4.io) 

Zn\ Zp Z n J 

In Figure [31 we used equation (|4.8I) to plot (f>* n as a function of cj>* . This curve is seen 
to describe very well the nonlinear saturation of the injected mass with increasing 
sample concentration in the reservoir. 

The denominator in equation (|2 . 15|) vanishes if — 9 /a. Far from the inlet, 
the sin gularity is approached a s — > 1/a (if a > 0), since, 9 ~ 1. It was shown 



earlier (jGhosal fc Chenl . l2010al) that the requirement of positivity of all ionic con- 

(4.11) 



centrations imposes a condition of self-consistency, <fi(x, t) < C where 
(z p - z„)/ {z p - z) if z > 



-[z n (z p - z n )]/[z p {z - z n )] if z < 0. 

Since, <p c < 1/a, if a > (see Appendix), the singularity is never reached if ini- 
tially 4> is less than (j> c everywhere in the channel. This, however, leaves open the 
possibility that the initial condition 0(0, t) might be such, that, exceeds C in 
some parts of the domain. We have now shown that if electrokinetic injection is 
used to insert the sample into the capillary, the shock at the inlet ensures that the 
injected sample concentration does not exceed §* m < 4>c (see Appendix), no matter 
how high the sample concentration in the reservoir may be. In capillary zone elec- 
trophoresis, one chooses the peak shape after electrokinetic sample injection as the 
"initial condition". Then, the only initial conditions possible are the "realizable" 
kind, where, 0(x, 0) < C . 

The speed of propagation of the analyte front in Figure [4] may be readily 
calculated when the front is far from the inlet. To do this, we note that since 
Pe = 200 ^> 1 and Ml, equation (|2.22[) may be written as 

where 

= 7~~~r ( 4 - 13 ) 

1 — am 



Then, the shock propagation speed is given by (IWhithaml . [Hzl the Rankine- 
Hugoniot jump condition 

Va = Q(M-Q(<h) (4 . 14) 
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where the subscripts 2 and 1 refer to conditions just behind the shock and just ahead 
of the shock respectively. Putting <fii = and <p2 = 4> m: we find, for the conditions 
shown in the lower panel of Figure El <$>* m = 3/7 « 0.43, so that, = 1.405 v. 
The shock speed determined from direct measurement of the front displacement in 
Figure [4] is 1.4 v, in close agreement with the theoretical prediction. In the upper 
panel, (j>* = 0.01, so that <t>* m « (j>* . In this case, both the theoretical and the 
measured values give V s ~ v . 



5. Conclusions 

The problem of electromigration of a sample ion from a large reservoir into a con- 
fined channel in the presence of a fully dissociated background electrolyte was stud- 
ied in the idealized situation where all ionic species have the same diffusivity. If the 
sample ion concentration is small compared to that of the background electrolyte, 
the behavior of the system is exactly as one might expect. When the electric field 
is switched on, a sample plug with concentration equal to that in the reservoir, is 
drawn into the microchannel. However, if the sample ion concentration is compa- 
rable or large compared to the background electrolyte, the behavior of the system 
is dominated by the intrinsic nonlincarity of the electro-diffusional problem, and 
the outcome is quite counter-intuitive. Once again, the sample is drawn into the 
microchannel as a plug of increasing length, but the concentration of sample ions in 
the plug is less than in the reservoir. In fact, as the reservoir concentration is pro- 
gressively increased, the concentration of sample ions in the capillary approaches a 
limit that depends on the valences of the three ionic species but is independent of 

the reservoir concentra tion. 

In an earlier paper dGhosal fc Chenl . l2010alV the authors presented a theory of 



electromigration of a sample in an infinitely long micro-channel when only three 
ions of equal diffusivity are present. There, it was shown, that, the requirement that 
none of the three ionic concentrations could ever be negative, put a limit to the 
validity of the theory. It was therefore required that the initial condition be such, 
that, the dimensionless sample concentration not exceed a certain (positive) critical 
concentration, (f> c . Here we have shown that initial conditions that do not satisfy 
this requirement cannot arise if the analyte is introduced into the micro-capillary 
by elecrokinetic injection. In fact, the diffusional boundary layer at the channel 
entrance ensures that the concentration of analyte drawn into the micro-capillary 
is always less than 4> c . Thus, the critical concentration can never be exceeed if the 
sample is introduced by electrokinetic injection, and consequently, the singularity 
at <f) = 1/a (if a > 0), inherent in equation (I2.22[) is never reached. 

One may question whether the strongly nonlinear regime considered here is 
of relevance to actual laboratory practice. The answer depends on the numerical 
value of the critical concentrations, and <p c . If the sample and carrier ions have 
similar valences, then these critical concentrations are all of order unity. Thus, 
to approach these critical values, the sample ions in the injected plug will need 
to be present at concentrations approaching that of the carrier electrolyte. Such 
high concentrations are normally not employed in laboratory practice with capil- 
lary electrophoresis. However, if the sample is a macro-ion the critical values may 
actually be qu ite small. For example , at pH 2.0 Bovine serum albumin has a va- 



lence, z ~ 55 (jFord fc Winzorl L1982). Then, in a univalent carrier electrolyte we 
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have 4> c ~ 0.04, so that the strongly nonlinear regime studied here may be easily 
reached. 

The shock like transition in concentration at the entrance of the micro-channel 

observed here is reminiscent of "de-sali nation shocks" that propagate outward from 

micro -channel / nano-channel junctions ( Mani. Zangle fc Santiagoll2009t IZangle. Mani fc Santiagol 
2009) or from the surface of a perm-selective membrane embedded in a microflu- 



idic channel. The mechanism of these shocks is related to the nonlinear coupling 
between the ion depletion zone (due to concentration polarization) near the mem- 
brane with the electric Debye layers at the walls of the micro-channel. The problem 
considered here is much simpler; Debye layers or concentration polarization are 
not involved. Nevertheless, it is another example of shock-like behavior that may 
be traced to the nonlinear nature of the underlying Nernst-Planck-Poisson sys- 
tem of equations for ionic transport. As in the case of de-salination shocks, the 
effect d escribed here would als o be relevant for solutes electromigrating into porous 
media ([Mani fc Bazanu 1201 If) or electromigration of analytes through var iable ge- 
omet ry channels, such as pores in membranes, where lubrication theory (|Ghosall . 



20021 ) can be used to generalize the analysis presented here. Finally, if the channel 



walls are charged, electroosmotic flow would arise so that the transport problem 
becomes intrinsically two or three dimensional. In such cases, homogenization can 
be used to achieve a reduced description, as shown by iGhosal fc Cher] (|2012h for 
electromigration in an infinitely long uniform channel. 
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Appendix: Proof of the inequality (p* m < <p c < 1/ a 
If a > 0, z p > z > z n . Then 

a<t> c =l Zp ~ z / n i Zp 7 T Zn l Mt ~, Zn , 7 (A-1 

\ -^-p ^ ^ TL ( %p ~ %Tl ) Zp 

thus, (f> c < 1/a if a > 0. 

To prove the first part of the inequality, (f>* n < <j) c , we first show that r > —z/z n . 
This is clearly true if z < 0. To show this when z > 0, first we use the electro- 
neutrality condition to express <j)* in terms of the other variables: 

d>* P = --r--K = -(-z-rz n ). (A-2) 

Zp Zp Zp 



Article submitted to Royal Society 



The nonlinear electromigration of analytes into confined spaces 



15 



Now we must have 0* > 0. This is always true if z < 0, but if z > 0, we require 
that r > —z/z n . 

We will now show that §* m < <p c by considering the two cases z > and z < 
separately. First suppose that z < 0. Then 

' < 7 A? \ = ZJL ^ = <t>c (A-3) 



m r+(z p - z)/(z p - z n ) (z p - z) I (z p - z n ) 
Now suppose that z > 0. Then 

1 



r 



+ ( Z P ~ Z )/( Z P - z n) 



(A-4) 



-(z/z n ) + (zp - z)/(z p - z n ) Zp z z n 
Thus, in all cases, <j>* < 4> c which completes the proof. 
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